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'-£^ ' Abstract: Using the spherically symmetric inhomogeneous Lemaitre-Tolman-Bondi dust 

solution, we study how the shear and the backreaction depend on the sharpness of the 
spatial transition between voids and walls and on the size of the voids. The voids considered 
here are regions with matter density Qq ~ and expansion rate H^to ~ 1, while the walls 
are regions with matter density Qq — 1 an d expansion rate -ffo^o — 2/3. The results 
indicate that both the volume-average shear and the variance of the expansion rate grow 
proportional to the sharpness of the transition and diverge in the limit of a step function, 
but, for realistic-sized voids, are virtually independent of the size of the void. However, the 
backreaction, given by the difference of the variance and the shear, has a finite value in the 
step- function limit. By comparing the exact result for the backreaction to the case where 
the shear is neglected by treating the voids and walls as separate Friedmann-Robertson- 
Walker models, we find that the shear suppresses the backreaction by a factor of (rn/io) 2 > 
the squared ratio of the void size to the horizon size. This exemplifies the importance of 
using the exact solution for the interface between the regions of different expansion rates 
and densities. The suppression is justified to hold also for a network of compensated voids, 
but may not hold if the universe is dominated by uncompensated voids. 
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1. Introduction 

An unresolved issue in cosmology is the effect of the structure formation on the cosmological 
observations beyond perturbative analysis of the Friedmann-Robertson- Walker or FRW 
models [1-4]. On account of the increased precision of these observations and the cosmic 
acceleration that they seem to indicate [5-7], the evaluation of this effect has become 
important [8-11] . A way to estimate the effect of the cosmic structures is via a backreaction 
term that arises by averaging inhomogeneous scalar quantities on spatial hypersurfaces [12]. 
The backreaction thus obtained is given by the variance of the expansion rate minus the 
non-negative average shear. 

A simplification in some estimates of the cosmological backreaction is based on partially 
or fully neglecting the shear on the interface between regions of different expansion rates 
[13-15]: these studies have found a significant amount of backreaction. On the other 
hand, perturbative studies that do not neglect the shear have suggested the backreaction 
to be insignificant [16-19]; however, see Ref. [4] for a recent discussion on the possible 
shortcomings of the perturbative approach in the backreaction problem. In this paper, we 
consider the issue within exact general relativity by studying: 



1. How the sharpness of the spatial transition between the voids (where ^o — and 
i^oio — 1) an d the walls (where £\) ~ 1 and Hoto ~ 2/3), as well as the size of the 
voids, affect the cosmic shear and thus the backreaction in the spherically symmetric 
inhomogeneous Lemaitre-Tolman-Bondi or LTB dust models. 

2. How the exact backreaction of the void- wall configuration relates to the case where 
the shear is neglected by treating the voids and the walls as separate FRW solutions. 

We perform the calculations on a spatial hypersurface at t = to and do not explicitly 
consider time-dependence (though to can be considered arbitrary). 

The systematic study on the role of shear makes our approach different from the 
previous studies on the dynamical backreaction in the LTB model which have focused on 
finding profiles that exhibit acceleration of the average expansion [20, 21] or on general 
properties of the backreaction [22, 23]. 

The paper is organized as follows. The necessary background of the LTB solution 
and the Buchert averaging method are introduced in Sects. ^ and ||, respectively. In Sect. 
§, we consider the shear and, in Sect. ||, the backreaction for almost compensated LTB 
voids residing within almost FRW walls. To show that the results do not depend on the 
compensating overdensity, we discuss uncompensated voids with a monotonically increasing 
physical matter density profile in Sect. |5.4| . Finally, the conclusions are given in Sect. |j 

2. LTB solution 

The exact spherically symmetric dust solution of general relativity was discovered by 
Lemaitre in 1933 [24] and is now commonly known as the LTB metric: 

ds 2 = _ df 2 + [^M)] 2 dr 2 + A 2 ( t){dQ 2 + ^2 M 2) (21) 

1 — k{r) 

where A'(r, t) = d r A(r, t), k{r) is related to the curvature of the spatial sections and A(r, t) 
is determined by the Friedmann-like evolution equation [25] 



H(r,t)=H (r) 
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A(r,t)J v uy " \A(r,t) 



(2-2) 



where H(r,t) = d t A(r,t)/A(r,t) = A(r,t)/A(r,t), H (r) = H(r,t ) and Q (r) are bound- 
ary condition functions specified on a spatial hypersurface t = to that determine the radial 
inhomogeneity profile, while the freedom to choose the function Aq{t) = A(r, to) corre- 
sponds to the scaling of the r-coordinate, used in this work to set 

A (r) = r . (2.3) 

The curvature function k(r) in the metric (|2.l| ) is related to these by 

k(r)^Hi(r)(no(r)-l)A 2 (r), (2.4) 



and the boundary condition function f2o( r ) is related to the physical matter density p(r,t 
on the t = to hypersurface as 

8ttG / B ( r ) Po(r)d 3 x 



n (r) 



3H 2 (r) f BM dh 



(2.5) 



where po(r) = p(r,to). Inversely, po{r) can be written in terms of fio( r ) anci Ho{r) as 



Po(r) 
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We only consider LTB models where the boundary condition functions obey the con 
straint 



1 



t y 3 

which approximates the simultaneous Big Bang condition 



H (r) = i 



v/i-^oM - tt (r)arsinh^kg« 
(l-^o(r)) 3 / 2 



(2.7) 



(2.8) 



such that the |error| < 1.5% in the considered interval < f^o < 1 and no error at the 
extremes Oo = and Oo = 1. We use Eq. ( |2.7D instead of Eq. ( |2.8D to make analytic cal- 
culations possible. The models with (approximately) simultaneous Big Bang form perhaps 
the most relevant subcase of LTB solutions, because, in these models, the inhomogeneities 
are growing modes (see [26]) as e.g. the near isotropy of the CMB suggests is also the case 
in the real universe. 

To study the shear in Sect. f| and the backreaction in Sect. ||, we need the following 
quantities of the LTB model: the shear scalar 

Jr. * - ^„„ = I f A (r,t) _ A'(r,t) N 2 



the volume expansion scalar 
6{r,t 



a r a 
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A(r,t) A'(r,t) 



(2.9) 
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A(r,t) A'(r,t) A'(r,t)A 2 (r,t) dr 
and their expressions on the t = to hypersurface: 

a 2 (r,to) = \{rH'o(r)) 2 , 



A 2 (r,t)A(r,t) 



(r,to) = 3H (r) + rH (r) 



l_d_ 



J ^o(r)) • 



Using the relation ( p.7| ), the shear (J27TTJ) can also be expressed in terms of Qo(r) as 



<j 2 (r,t ) 



1 (rV (n (r)) 2 
54 \to) fio(r) 
We also have use for the LTB volume element 

/gdrdedy=^£M^-drd9d<p 

where g denotes the determinant of the spatial part of the metric. 



(2.10) 

(2.11) 
(2.12) 

(2.13) 

(2.14) 



2.1 The void- wall model 

We use the LTB solution to model a configuration that consists of two different regions: a 
void where Qq ~ and a wall where f^o — L with a smooth transition in between. For this 
we choose the £lo( r ) profile as follows: 

n (r) = (l-e-^ n ) 2 , (2.15) 



where r^ determines the size of the void and n the sharpness of the spatial transition 
between the void and the wall. By virtue of the simultaneous Big Bang condition ([LTj), 



Eq. (2.15) implies 



ff (r) = ^(2 + e-OV*O n ) , (2.16) 

which tells us that the void expands faster than the wall by a factor of 3/2. 

Using LTB models with the void-wall profile (2.15), we can easily join together many 



voids to construct a model for a network of voids. This is because LTB solutions with 
the profile fl2,15| ) are, up to negligible terms of order e~( R ' r °> , $7 = 1 FRW dust solutions 
outside the void at R > tq, all with to as the age of the universe, implying that the different 
LTB solutions naturally join together in the wall region. 

2.2 Small k(r) expansion 

For a realistic or sub-horizon size void r^ <C to, which we next show implies \k(r)\ <C 1 V r. 
Note, however, that since the Ricci scalar of the spatial sections is given by 

^'(r.oV*) !^''*''""- (2 ' 17) 

or at t = to by 

no = ^^(rk(r)), (2.18) 

the condition \k(r)\ <C 1 does not imply that the spatial curvature is small in the void. 



With the conditions (2.3) and (2.7), the curvature function ( |2.4p becomes 

k ^ = -{£) (i-^V^M) (1-Oo(r)). (2.19) 

* v ' 

<1 

For the considered profiles, £lo(r) rapidly approaches the value f^o = 1 outside the void or 
when r > ro, thus implying 

max(|/c(r)|) ~ ( — J <C 1 . (2.20) 



As is evident in Eqs. (|2.4| ) and ( 2.19|) , k{r) depends only on the inhomogeneity profile 



functions fio(r) and/or Ho(r), but not on their derivatives. Hence the small k(r) approx- 
imation is equally valid for all the profiles considered here regardless of the sharpness of 
the transition between the void and the wall. 



Since \k(r)\ <C 1, the part of the LTB volume element ( 2.14 ) that contains the curvature 
function k(r) can be expanded as follows 

l + \k{r) + 0{k 2 {r)) . (2.21) 



VT^W) 2 

Furthermore, Eqs. ( |2.ip , ( |2.3| ) and fl2.20| ) imply that the coordinate r closely measures the 
proper distance on the spatial hypersurface defined by t = t . 



3. Scalar averaging 

The spatial volume-average of a scalar S(x,t) is defined as 



f v S(x, t)y/gd 



,r 



{s(x,t)) v = j » : '/^ — , (3.i) 

Jv V9 a x 

where g is the determinant of the spatial metric and T> is the averaging domain or a region 
of the spatial sections. 

By applying the averaging ( |3.1| ) to the scalar parts of the Einstein equation in an 
irrotational dust universe, we obtain the Buchert equations [12]: 

3^ = -*xG{p)„(t) + Q v (t) (3.2) 

av(t) 

3 (SH) =^ G ^{t)-\{K)v(t)- l -Q v {t) (3.3) 

dt{ P )v(t) = -3^§-(p)v(t) , (3.4) 

where the regional scale factor is defined as 

^ ( j v ^g{x\t)tfx \ V3 
\/x>V / 5(^Vo)d 3 a;y 

and the backreaction term 

Gd(*) = f ((^ 2 )7. " («>?,) - (^)t. (3.6) 

quantifies the difference of the time evolution of the averages relative to the homogeneous 
quantities in an FRW dust universe. In what follows, we take the averaging domain to be 
an origin-centered ball of coordinate radius R, i.e. V = M(R), and do not write it explicitly 
anymore, so from here on for all scalar s S: 

(S)v = (S) M(R) = (S) . (3.7) 



3.1 Averaging in the LTB model 



Given the LTB volume element ( |2.14 ), we can expand Eq. ( |3.1[ ) for small k(r) using the 
result ( ggg ): 

(3.8) 



(S) = (5) + I ({Sk) - (S) (k)o) + 0(k 2 ) , 



where the subscript now refers to averages where k(r) = in the integration measure, 
that is y/go — ^4'( r ) t)A 2 (r, t) sin 8. Using this notation, the expression of Qo( r ) m Eq. ( |2.5| ) 
can be recognized simply as 

<Po(r))o 



n n (r) 



Pcrit(r) 



(3.9) 



where p cvit (r) = 3# 2 (r)/8vrG. 



4. Shear 



When the transition between the void and the wall is made sharper by increasing n, the 
shear density on the interface grows but it also becomes more localized. This can be seen in 
Eq. ( 2.13| ) and is illustrated in Fig. [I], where we plot a 2 (r, to) for a few profiles ( [2.15 ) with 
different values of n. It is hence a priori unclear how the integrated or volume-averaged 
shear behaves when n is varied. 
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Figure 1: Shear distribution as a function of r for three different values of n: n = 2 (red dash dot 
curve), n — 4 (black dashed curve) and n = 6 (blue solid curve). All profiles in this figure have 
r = 0.01to. 



We apply the expansion (|3.8|) for the shear: 

<^ 2 > = <^ 2 >o + \ ((a 2 k} - (a 2 ) (k) ) 



(4.1) 



=(^ 2 h 



Since the shear falls off rapidly in the wall outside the void, we can also make the following 
approximation 



R roo 

a (r,to)r dr ~ / a (r,to)r dr 

n Jo 



(4.2) 



where R > r$ is the (coordinate) radius of the spherical averaging region. The validity of 
the approximation fl4.2|) depends on n: n = 1 requires R/vq > 6, n = 2 requires R/tq > 2, 
n = 6 requires RJTq ^1.2, while for very large n it is enough just to have R/tq > 1; see 
Fig. | 
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Figure 2: The area above the curve (conservatively) represents the region in the parameter space 
where the approximation (jh2) is valid. 



Using Eqs. ( 2.13| ), Q2.14 ), ( 2.15J ) and (3.1), along with the approximation (|4.2j), we 
obtain the following analytic expression for the average shear in the zeroth order of k(r): 



<^(§) 3 R) r (!'-" 



(4.3) 



where T stands for Euler's gamma function. Similarly, the first order term for shear is 
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n / V Ti 
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r [ - ) 2~ 3 / n [-8 - 4 • 2~ 5 / n + 2 • 3" 5/n + i~ 5/n 



which is suppressed by the overall factor (ro/io) 2 relative to the leading order term (|4.3|). 
Therefore, as can be easily verified by numerical methods, for the values of R/vq and n 
above the curve in Fig. |2[ Eq. (|4.3| ) gives the correct result for the average shear. 

An immediate conclusion from the expression ( |4.3D is that, for realistic or sub-horizon 
size voids, the average shear is independent of the size of the void (as long as the volume- 
ratio (ro/i?) 3 is kept fixed for each void- wall pair). Since, as explained in Sect. 2.1, we can 



join together many LTB void-wall profiles fl2.15| ) to form a network of voids, Eq. (4.3) thus 
also gives the average shear for the network of sub-horizon size voids where each void can 
have a different value of r$ . 

For a universe with half of the volume in voids 1 , (tq/R) 3 = 1/2. Using this value and 
setting n = 5 in Eq. (|4.3|), we obtain 



In the step function limit n 



<^ 2 Mo)> 

► oo, Eq. 



O.ltr 



(4.5) 



<-V,to)>~^ - 2 (?)V 



1 

18~ 



diverges with the following asymptotics: 

(4.6) 



RJ 



telling us that the sharper the transition between the void and the wall, the higher the 
value of the average shear. 

5. Backreaction 



(4.4) 



We consider here the backreaction (3.6) for LTB models with a void-wall profile given by 
Eqs. ( [OBD and (|2l6| ). 



is helpful to notice that the following terms can be written as total derivatives: 



5.1 Backreaction in Oth order of k(r) 

Let us first show that the backreaction (p.6|) vanishes in the Oth order of k{r). For this it 



(5.1) 
(5.2) 







3 ' ~ A 2 AA> ~ A 2 A'dr 



A 2 A 



A A' _ 1 d 
A + A J ~ A 2 A'~ifr 



AA' 



Since the volume element ( p. 14 ) in the Oth order of k{r) is just 



1 2 /i / 



-gQ~& A x = AM'sin0drd0dp , 



(5.3) 



1 Order of magnitude consistent with observed values quoted e.g. in Ref. [27]. 



we have 



and 



so 



^-^-isi^efjlf (5.4) 

Co = J «« 2 >o - WS) - {°*h = (|» 2 - »*)„ - \{»)l = • (5.6) 

Given the zeroth order term of the backreaction ( |5,6| ) vanishes, the variance of the 
expansion rate must be equal to the average shear ( |4.3| ) to the leading order, so we have 

£ (< ^ _ WS) = £ (£)* (x + 5) r (1)**.. (5.7) 

5.2 Backreaction in 1st order of k(r) 

As the average shear and the variance of the expansion rate cancel exactly in the leading 
order of k(r), we must go beyond the zeroth order to obtain the leading order term of the 
backreaction. 

Using the result (J3iJ|) , the backreaction (|3.6| ) becomes 

"2 



Q = I ((9 2 )o - (9)1) - (A, + (k)o (\(9)l - \ 



§«# 2 >o - (0)1) - (a 2 ) 



+ 



+ ( (t " t) *>. - f <(,> ° <«*>• < (5 - 8) 

where, according to the result (fT6|), Qo = | ((# 2 )o — (#)o) ~~ (°" 2 )o vanishes, so Eq. ( p^ ) 
reduces to 

S = 3Hl(R)(k) + ((!L-^\ k) Q - 2H (R)(9k) , (5.9) 



where we have written (0)o in terms of Hq(R) using the result (|5.5| ). 

Using the approximation ( |4.2[) for the integrals in Eq. (|5.9| ) that contain fe(r) as a 
common factor, we obtain 

Q = V (-) 3 f-V -r f ^ (* • s- 5 /" - 1. 4- 5 /" + 1. 5- 5 /" + - • 6- 5 /"+ 

9° Vfl/ V^o/ n Wl 9 3 3 9 

+ e -2{R/r r C_8 - 4 ■ 2" 5/n + 2 • 3~ 5/n + 4~ 5/nN ) + 



_ e -(«/ro)» ( ° . 2 "5/n + 32 3 _ 5/n _ 7 4 _ 5/n _ 4 5 „ 5/n j (5 1Q) 

3 9 3 3 



In the limit n — > oo, Eq. ( 5.10J ) gives 



«-^®'(S)"- (5I1) 



4 ,-2 / r 0\ 3 f H) 

showing that although both the average shear and the variance of the expansion rate diverge 
as n — > oo, the backreaction, given by their difference, has a finite limit. Note that, in the 
validity region of Fig. y, the backreaction ( |5.10 ) is qualitatively similar for all values of 



n. 



but we consider the step function limit to make comparison with the model that consists 



of two disjoint FRW solutions in Sect. 5.3 



5.3 Comparison with a model that consists of two disjoint FRW solutions 

Let us compare the exact solution of the void- wall configuration to a model where the walls 
and voids are treated as separate FRW regions, following the idea outlined in Ref. [13]. 
This comparison directly addresses the role of shear since the shear is by construction zero 
in the disjoint FRW approximation. 

For the void or the £1 = region we use the Milne solution for which 

3 



t 



(5.12) 



and for the wall or the fi = 1 region we use the Einstein de Sitter or EdS solution for which 

2 



t 



(5.13) 



Since the shear vanishes in the FRW models, the backreaction is given just by the variance 



Qfrw = 3 {(0 2 ) ~ (0) 2 ) , 



where 



and 



if 



(0) 



v w + v v 

V W 6 W + V V 9 V 



v w + v v 

By inserting Eqs. ( |5.12| ) and ( 5.13[ ) in Eqs. ( |5.15| ) and ( |5.16 ), we obtain 

2 



Qfrw = -^t Q /v(l — /«) > 



where f v is the void fraction: 



fv 



V v 



v v + v w 



For realistic-sized voids, r$ <C to, we have 

4 



V v = gvrrg 






3 ° 



while the wall is spatially flat, so Eq. ( 5.17 ) can also be written as 

2 „ /rn\3 T /Tn\3 



«™-M3) -S) 



(5.14) 

(5.15) 
(5.16) 

(5.17) 
(5.18) 

(5.19) 
(5.20) 



For the void fraction value f v = 1/2 or R = 2 1 ' 3 ro, the backreaction (|5.20| ) gets its 
maximum value: 

(5.21) 



Qfrw — ^^o 



1 

-i 
6 



We can compare this to the LTB result by setting R = 2 ' 3 ro likewise in Eq. ( 5.11| ): 

2 



Q 



135 



-t. 



o 



2 ( '[O 

to 



(5.22) 
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implying the discrepancy: 



q = ^(t ] tJnm 



(5.23) 



The result ( |5.23| ) demonstrates how important it can be to take into account the shear: 
for a realistic-sized void with tq = 0.01£o> the FRW approximation overestimates the 
backreaction by the tremendous factor of 10 5 . The suppressive factor (ro/io) 2 appears 
to be consistent with the results from perturbative analysis in Refs. [16, 19]. 

5.4 Uncompensated voids 

The void- wall profiles ( 2.15 ) considered in Sects. I-5J3 have a compensating overdense 



peak in the matter density between the void and the wall regions; see Fig. |3|. The peak is 
required to make the wall (very close to) a flat FRW region, as can be seen from Eqs. fl2.£| ) 
and ([^Ij). In this section, we consider profiles without the peak to see whether it plays an 
important role in suppressing the backreaction relative to the FRW value as in Eq. ( 5,23; ). 
Instead of specifying &o(r) and to, a simultaneous Big Bang LTB model can alterna- 
tively be parameterized by the physical matter density profile po(r) on a spatial hyper- 
surface at t = to- In this way, we can avoid the compensating over density by choosing a 
monotonically increasing density profile such as 



MO = Pcri t (oc)(l-e-^)") 



(5.24) 



where /0 cr ; t (oo) = 1/QitGtQ. Substituting Eqs. ( f2~7^ ) and ( |5.24 ) in Eq. ( |2.5[) , we obtain 

2 



Sh(r) = - I I 



N 



1 



3V3 



E 



-1)' +1 (r/ro) 



111 



1 l\{nl + 3) 



(5.25) 



which, by virtue of Eq. (2.7), implies 

.1/ 
\ 



tr 



Ho(r) = ^ 



1 + 



N 



3V3 



E 



l\(nl + 3) 



(5.26) 



In the step function limit n — > 00, Eqs. ( 5.25J ) and ( ^.26 ) become 



fio(r) = ^ I J 



fWW)'«'-«> 



(5.27) 



and 



Ho(r) = to" 1 



k^-fWW)^-^ 



(5.28) 



where @(r — tq) = lim n _ >00 (l — e ' r / r °) n ) stands for the Heaviside step function. 
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The backreaction ( |3,6| ) cannot be integrated analytically for the profile ( |5.24 ). Instead, 
performing the integrals numerically for the parameter values (ro/R) s = 1/2 and vq = 
0. Olio , we obtain 

Q = 1.6 • KT 6 ^ 2 . (5.29) 



This differs very little from the case with the step function in ilo( r ) an d the compensating 
overdensity in the matter distribution: for comparison, substituting the values (ro/R) 3 = 
1/2 and r = 0.01i in Eq. (|5T2p yields 



Q = 1.5 • 10~ b t, 



-2 



(5.30) 



The near equality of the results (|5.29| ) and ( |5.30|) implies that the suppression factor ap- 
pearing in Eq. ( |5.23D is not due to the overdense peak between the wall and the void 
regions. Furthermore, keeping the void-wall volume ratio fixed to (ro/R) 3 = 1/2 and vary- 
ing ro shows that the dependence of the backreaction on rg is very similar to Eq. ( 5.22J ), 



i.e. Q oc 7q also for the profile ( 5.24 ) 



6 Jt G t Q p Q 




0.03 



Figure 3: Matter density as a function of r for the profile (5.24) with n — > oo (magenta dotted 
curve) and the profile ( [2.15| ) with three different values of n: n = 2 (red dash dotted curve), n = A 
(black dashed curve), n = 6 (blue solid curve). All profiles in this figure have tq — O.Olto. 



Extrapolating the result from a void-wall pair to the network of voids is harder for 
the uncompensated voids than for the compensated ones. The reason is that, without 
the compensating overdensity, the solution approaches FRW metric slower: although the 
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matter density is constant outside the void (r > ro) for the profile ( 5.24 ) in the limit 
n — > oo, the spatial curvature becomes (nearly) constant only much further (r > 10ro) 
from the void. It thus appears that in order to match together uncompensated voids in 



the simple fashion explained in Sect. 2.1, the separation between the voids must be so 
large that only a few voids fit inside the horizon. Therefore, more sophisticated junction 
conditions need to be applied in joining together uncompensated voids to obtain a global 
void- fraction more consistent with observations such as the ones quoted in Ref. [27]. 



6. Conclusions 

We have considered the role of shear in the cosmological backreaction problem using the 
spherically symmetric LTB dust solution. The LTB models utilized in Sects. ^-5^ have 



a close to simultaneous Big Bang with the inhomogeneity profile (2.15) that at t = to 
interpolates between a void region (where ^o — and H^to ~ 1) and a wall region (where 
Qq ~ 1 and Hoto ~ 2/3) with the parameters ro and n describing the size of the void and 
the sharpness of the transition between the two regions, respectively. 

For realistic-sized voids (ro <C to) with a fixed void-wall volume ratio, the results of 
Sect. |] show that the volume- average shear is independent of ro or the size of the void. 
This, along with the fact that the different void- wall solutions naturally join together in 
the EdS-like wall region, implies that the volume- average shear found for a single void-wall 
pair actually applies also for a network of voids of different size. Moreover, we found that 
for values n > 10 both the volume-average shear and the variance of the expansion rate 
grow proportional to n, hence diverging in the limit n — > oo. However, as shown in Sect. [5], 
the backreaction, given by the difference of the variance and the shear, has a finite limiting 
value for n — > oo. 

In Sect. |5.3j , to estimate the role of the shear in the backreaction, we compared the 
LTB result for the backreaction in the limit n — > oo to a simplified model that ignores 
the shear by approximating the void- wall configuration with two disjoint FRW solutions: 
Slo = or the Milne solution for the void and f&o = 1 or the EdS solution for the wall. 
The comparison shows that the backreaction obtained using the exact LTB solution is 
suppressed by at least a factor of (ro/to) 2 relative to the value obtained from the model 
made up of the separate FRW solutions, hence implying that the shear plays a crucial role 



here. In Sect. 5^, we demonstrated that the suppression is not due to a compensating 
overdensity between the void and the wall regions. 

The LTB-based void-wall models considered here contain various simplifications that 
may misestimate the backreaction of the real universe. Firstly, the models lack globally 
significant spatial curvature which is known to be important for large backreaction [9]. 
Large negative spatial curvature can arise in the late universe if its volume becomes domi- 
nated by uncompensated voids. This can be seen by realizing that a dense enough network 
of uncompensated voids is dynamically similar to a huge negatively curved void of size 
r oAo — C(l)- Another issue is the spherical symmetry of the voids which is obviously 
broken in the real universe. Finally, the models considered here do not contain collapsing 
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regions, which have been suggested to play an essential role in the average dynamics of the 
universe [13, 14]. 

Besides the dynamical backreaction ([O]) studied in this work, the cosmological struc- 
tures may have other effects on the observations that are not captured by a spatial averaging 
procedure. Proposed examples include effects on the propagation of light [10, 28-32] and 
effects due to our non-average location in the universe [11,15,33,34]. Therefore, small 
dynamical backreaction does not alone imply that the effects of the inhomogeneities on the 
observations were small or insignificant. 
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